Search for: All records

Abstract At the edge of alpine and Arctic ecosystems all over the world, a transition zone exists beyond which it is either infeasible or unfavorable for trees to exist, colloquially identified as the treeline. We explore the possibility of a thermodynamic basis behind this demarcation in vegetation by considering ecosystems as open systems driven by thermodynamic advantage—defined by vegetation’s ability to dissipate heat from the earth’s surface to the air above the canopy. To deduce whether forests would be more thermodynamically advantageous than existing ecosystems beyond treelines, we construct and examine counterfactual scenarios in which trees exist beyond a treeline instead of the existing alpine meadow or Arctic tundra. Meteorological data from the Italian Alps, United States Rocky Mountains, and Western Canadian Taiga-Tundra are used as forcing for model computation of ecosystem work and temperature gradients at sites on both sides of each treeline with and without trees. Model results indicate that the alpine sites do not support trees beyond the treeline, as their presence would result in excessive CO$$_2$$ ${}_2$loss and extended periods of snowpack due to temperature inversions (i.e., positive temperature gradient from the earth surface to the atmosphere). Further, both Arctic and alpine sites exhibit negative work resulting in positive feedback between vegetation heat dissipation and temperature gradient, thereby extending the duration of temperature inversions. These conditions demonstrate thermodynamic infeasibility associated with the counterfactual scenario of trees existing beyond a treeline. Thus, we conclude that, in addition to resource constraints, a treeline is an outcome of an ecosystem’s ability to self-organize towards the most advantageous vegetation structure facilitated by thermodynamic feasibility. 

Abstract. Wetlands are the largest natural source of methane (CH4) emissions globally. Northern wetlands (>45° N), accounting for 42 % of global wetland area, are increasingly vulnerable to carbon loss, especially as CH4 emissions may accelerate under intensified high-latitude warming. However, the magnitude and spatial patterns of high-latitude CH4 emissions remain relatively uncertain. Here, we present estimates of daily CH4 fluxes obtained using a new machine learning-based wetland CH4 upscaling framework (WetCH4) that combines the most complete database of eddy-covariance (EC) observations available to date with satellite remote-sensing-informed observations of environmental conditions at 10 km resolution. The most important predictor variables included near-surface soil temperatures (top 40 cm), vegetation spectral reflectance, and soil moisture. Our results, modeled from 138 site years across 26 sites, had relatively strong predictive skill, with a mean R2 of 0.51 and 0.70 and a mean absolute error (MAE) of 30 and 27 nmol m−2 s−1 for daily and monthly fluxes, respectively. Based on the model results, we estimated an annual average of 22.8±2.4 Tg CH4 yr−1 for the northern wetland region (2016–2022), and total budgets ranged from 15.7 to 51.6 Tg CH4 yr−1, depending on wetland map extents. Although 88 % of the estimated CH4 budget occurred during the May–October period, a considerable amount (2.6±0.3 Tg CH4) occurred during winter. Regionally, the Western Siberian wetlands accounted for a majority (51 %) of the interannual variation in domain CH4 emissions. Overall, our results provide valuable new high-spatiotemporal-resolution information on the wetland emissions in the high-latitude carbon cycle. However, many key uncertainties remain, including those driven by wetland extent maps and soil moisture products and the incomplete spatial and temporal representativeness in the existing CH4 flux database; e.g., only 23 % of the sites operate outside of summer months, and flux towers do not exist or are greatly limited in many wetland regions. These uncertainties will need to be addressed by the science community to remove the bottlenecks currently limiting progress in CH4 detection and monitoring. The dataset can be found at https://doi.org/10.5281/zenodo.10802153 (Ying et al., 2024). 

Abstract Climate change is rapidly altering composition, structure, and functioning of the boreal biome, across North America often broadly categorized into ecoregions. The resulting complex changes in different ecoregions present a challenge for efforts to accurately simulate carbon dioxide (CO₂) and energy exchanges between boreal forests and the atmosphere with terrestrial ecosystem models (TEMs). Eddy covariance measurements provide valuable information for evaluating the performance of TEMs and guiding their development. Here, we compiled a boreal forest model benchmarking dataset for North America by harmonizing eddy covariance and supporting measurements from eight black spruce (Picea mariana)-dominated, mature forest stands. The eight forest stands, located in six boreal ecoregions of North America, differ in stand characteristics, disturbance history, climate, permafrost conditions and soil properties. By compiling various data streams, the benchmarking dataset comprises data to parameterize, force, and evaluate TEMs. Specifically, it includes half-hourly, gap-filled meteorological forcing data, ancillary data essential for model parameterization, and half-hourly, gap-filled or partitioned component flux data on CO₂(net ecosystem production, gross primary production [GPP], and ecosystem respiration [ER]) and energy (latent [LE] and sensible heat [H]) and their daily aggregates screened based on half-hourly gap-filling quality criteria. We present a case study with the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) to: (1) demonstrate the utility of our dataset to benchmark TEMs and (2) provide guidance for model development and refinement. Model skill was evaluated using several statistical metrics and further examined through the flux responses to their environmental controls. Our results suggest that CLASSIC tended to overestimate GPP and ER among all stands. Model performance regarding the energy fluxes (i.e., LE and H) varied greatly among the stands and exhibited a moderate correlation with latitude. We identified strong relationships between simulated fluxes and their environmental controls except for H, thus highlighting current strengths and limitations of CLASSIC. 

Abstract Northern peatlands are a globally significant source of methane (CH₄), and emissions are projected to increase due to warming and permafrost loss. Understanding the microbial mechanisms behind patterns in CH₄production in peatlands will be key to predicting annual emissions changes, with stable carbon isotopes (δ¹³C‐CH₄) being a powerful tool for characterizing these drivers. Given that δ¹³C‐CH₄is used in top‐down atmospheric inversion models to partition sources, our ability to model CH₄production pathways and associated δ¹³C‐CH₄values is critical. We sought to characterize the role of environmental conditions, including hydrologic and vegetation patterns associated with permafrost thaw, on δ¹³C‐CH₄values from high‐latitude peatlands. We measured porewater and emitted CH₄stable isotopes, pH, and vegetation composition from five boreal‐Arctic peatlands. Porewater δ¹³C‐CH₄was strongly associated with peatland type, with δ¹³C enriched values obtained from more minerotrophic fens (−61.2 ± 9.1‰) compared to permafrost‐free bogs (−74.1 ± 9.4‰) and raised permafrost bogs (−81.6 ± 11.5‰). Variation in porewater δ¹³C‐CH₄was best explained by sedge cover, CH₄concentration, and the interactive effect of peatland type and pH (r² = 0.50,p < 0.001). Emitted δ¹³C‐CH₄varied greatly but was positively correlated with porewater δ¹³C‐CH₄. We calculated a mixed atmospheric δ¹³C‐CH₄value for northern peatlands of −65.3 ± 7‰ and show that this value is more sensitive to landscape drying than wetting under permafrost thaw scenarios. Our results suggest northern peatland δ¹³C‐CH₄values are likely to shift in the future which has important implications for source partitioning in atmospheric inversion models. 

Abstract Arctic wetlands are known methane (CH₄) emitters but recent studies suggest that the Arctic CH₄sink strength may be underestimated. Here we explore the capacity of well-drained Arctic soils to consume atmospheric CH₄using >40,000 hourly flux observations and spatially distributed flux measurements from 4 sites and 14 surface types. While consumption of atmospheric CH₄occurred at all sites at rates of 0.092 ± 0.011 mgCH₄ m⁻² h⁻¹(mean ± s.e.), CH₄uptake displayed distinct diel and seasonal patterns reflecting ecosystem respiration. Combining in situ flux data with laboratory investigations and a machine learning approach, we find biotic drivers to be highly important. Soil moisture outweighed temperature as an abiotic control and higher CH₄uptake was linked to increased availability of labile carbon. Our findings imply that soil drying and enhanced nutrient supply will promote CH₄uptake by Arctic soils, providing a negative feedback to global climate change. 

Abstract The Arctic–Boreal Zone is rapidly warming, impacting its large soil carbon stocks. Here we use a new compilation of terrestrial ecosystem CO₂fluxes, geospatial datasets and random forest models to show that although the Arctic–Boreal Zone was overall an increasing terrestrial CO₂sink from 2001 to 2020 (mean ± standard deviation in net ecosystem exchange, −548 ± 140 Tg C yr⁻¹; trend, −14 Tg C yr⁻¹;P < 0.001), more than 30% of the region was a net CO₂source. Tundra regions may have already started to function on average as CO₂sources, demonstrating a shift in carbon dynamics. When fire emissions are factored in, the increasing Arctic–Boreal Zone sink is no longer statistically significant (budget, −319 ± 140 Tg C yr⁻¹; trend, −9 Tg C yr⁻¹), and the permafrost region becomes CO₂neutral (budget, −24 ± 123 Tg C yr⁻¹; trend, −3 Tg C yr⁻¹), underscoring the importance of fire in this region. 

Abstract Process‐based land surface models are important tools for estimating global wetland methane (CH₄) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site‐level patterns of freshwater wetland CH₄fluxes (FCH₄) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model‐observation disagreements are mainly at multi‐day time scales (<15 days); (b) most of the models can capture the CH₄variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH₄production). Our evaluation suggests the need to accurately replicate FCH₄variability, especially at short time scales, in future wetland CH₄model developments. 

Environmental observation networks, such as AmeriFlux, are foundational for monitoring ecosystem response to climate change, management practices, and natural disturbances; however, their effectiveness depends on their representativeness for the regions or continents. We proposed an empirical, time series approach to quantify the similarity of ecosystem fluxes across AmeriFlux sites. We extracted the diel and seasonal characteristics (i.e., amplitudes, phases) from carbon dioxide, water vapor, energy, and momentum fluxes, which reflect the effects of climate, plant phenology, and ecophysiology on the observations, and explored the potential aggregations of AmeriFlux sites through hierarchical clustering. While net radiation and temperature showed latitudinal clustering as expected, flux variables revealed a more uneven clustering with many small (number of sites < 5), unique groups and a few large (> 100) to intermediate (15–70) groups, highlighting the significant ecological regulations of ecosystem fluxes. Many identified unique groups were from under-sampled ecoregions and biome types of the International Geosphere-Biosphere Programme (IGBP), with distinct flux dynamics compared to the rest of the network. At the finer spatial scale, local topography, disturbance, management, edaphic, and hydrological regimes further enlarge the difference in flux dynamics within the groups. Nonetheless, our clustering approach is a data-driven method to interpret the AmeriFlux network, informing future cross-site syntheses, upscaling, and model-data benchmarking research. Finally, we highlighted the unique and underrepresented sites in the AmeriFlux network, which were found mainly in Hawaii and Latin America, mountains, and at under- sampled IGBP types (e.g., urban, open water), motivating the incorporation of new/unregistered sites from these groups.